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Abstract 

We present a simple model for the dimuon production from the decay of vector 
mesons which are hadronized from a baryon-free quark-gluon plasma. Contin- 
uous decay of vector mesons in the medium is known to enhance the dimuon 
production and in our model we can estimate the total yield of dimuons de- 
cayed from vector mesons through the rate equations for the number of vector 
mesons of interest. Muon pair production critically depends on the time du- 
ration of both the mixed and hadronic phases. Results are discussed in view 
of possible measurements in experiments. 
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I. INTRODUCTION 



Recently dilepton production from vector meson decays in relativistic heavy-ion collisions 
has been drawing stong interest |l] . One of the reasons is that the properties of vector mesons 
such as masses or widths may change in a hot and dense medium as a result of the onset 
of chiral symmetry restoration |2| and dileptons under vector meson channels may serve as 
good probes of such changes ||,[| . 

The other reason for the interest is due to the fact that vector mesons continously 
decay in the hot and dense medium and thus they may contain valuable informations on 
the spacetime evolution of the collision zone. Heinz and Lee |5| suggested to measure the 
number of dileptons under rho and omega meson peaks from which one can estimate the 
lifetime of the purely hadronic matter. Seibert et al. suggested to measure the slope of the 
transverse mass spectrum of dileptons as a thermometer of the QCD transition ||. 

Or combining the obove two mechanisms there may appear double dilepton peaks for 
phi or rho mesons as a consequence of the expected quark-gluon plasma phase transition. 
However, in studying the dileptons one relies only on the idea and still important quantities 
such as relative magnitude of the two peaks are not estimable. 

In this paper we extend the model by Heinz and Lee || for the dilepton production from 
the purely hadronic matter to the case of the evolution of quark-gluon plasma and study the 
total yield of dileptons from the decay of vector mesons. This model was initailly used to 
simultaneously explain the J/\l/ supression and enhancement of <p/(p+uj) measured by NA35 
from a purely hadronic matter and modified to include the free decay to study the muon 
pair production ||. We believe it is a simple model for the number of resonant dileptons 
which allows easy modification to incorporate the change of vector meson properties in the 
medium. 

For the dilepton production decayed from vector mesons, V, we need to know the time 
evolution of the number of vector mesons, V, which can be obtained assuming a model both 
for the dynamic evolution such as hydrodynamics and for the hadronization mechanism. 
Hydrodynamics becomes simple for a baryon-free quark-gluon plasma whose physical quan- 
tities are invariant under longitudinal Lorentz boosts [7| and for the hadronization scenario 
we assume that hadronization occurs as a first order phase transition maintaining chemical 
equilibrium. These assumptions lead to a constant rate of hadronization. By assuming 
model equations of state one can evaluate the number of vector mesons as a function of time 
starting from an initial conditions reached during the collisions at a given beam energy. 

However, continuous decays of vector meson V may make deviations from chemical 
equilibrium and in this paper, we construct rate equations for the time evolution of vector 
mesons, V, by considering the gain term due creation of V via secondary collisions and loss 
terms due to annihilation through secondary collisions and free decay of V together with 
the hadronization rate in the mixed phase. One can solve the rate equations for the particle 
species of interest assuming that other particle densities are determined from the thermal 
and chemical equilibrium. In this way we loose knowledge on the time scales of the evolution 
and will present our results as a function of the time durations of each phases by controlling 
the degree of equilibration, namely the ratio of gain and loss terms. 

In section 2, the equilibrium hadronization of a Lorentz-boost invariant quark-gluon 
plasma is discussed. Hadronization at constant temperature implies a constant rate of 
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hadronization during the mixed phase. In section 3, considering secondary collisions among 
hadrons and the free decay of certain particle species, rate equations for the change in 
particle numbers are constructed and solved in both the mixed and hadronic phases. Muon 
production decayed from the particle of interest are formally obtained. In section 4, the 
results of our numerical calculations are presented and we show that dilepton production 
from rho meson decays depends critically on the time length of both the mixed and hadronic 
phase. We summarize in section 5. 



II. EQUILIBRIUM HADRONIZATION OF A BARYON-FREE QUARK-GLUON 

PLASMA 

In this section we discuss the hadronization rates of particles under the assumption 
that thermal and chemical equilibrium is maintained throughout the phase conversion of a 
baryon-free quark-gluon plasma. This assumption considerably simplifies the calculations, 
but our final argument does not qualitatively depend on it. The results depends crucially, 
however, on our assumption that the phase transition is of first order. This issue is still 
under debate [I]. 

The most important property of a baryon-free quark-gluon plasma which is created in 
the central region during relativistic heavy-ion collisions is the invariance of the rapidity 
distribution under longitudinal Lorentz boosts. This property makes the hydrodynamic 
description very simple @: 

d A S d A N 

— = and — — = 0. 1) 

dx A dx A v ; 

Assuming only longitudinal expansion (which initially dominates in any case [ref]) it is 
convenient to use the rapidity y and the longitudinal proper time r as varaibles. Then 
above hydrodynamic equations for the perfect fluid implies 

s( t )t = s(r )r , p(r)r = p(r )r . (2) 

Here p can be the number density of any particle species in the system, e.g. the quark 
density or the baryon density, etc. 

In other words, the spatial volume increases linearly in r as a consequence of hydrody- 
namic expansion: 

V(r) = -V(r ). (3) 

TO 

We assume Eq. (|2|) holds throughout the entire evolution of a quark-gluon plasma, namely 
through the mixed and hadronic phases until freeze-out. 

As the system cools below due to its hydrodynamic expansion the confinement phase 
transition occurs. We assume that this phase transition is of first order. Conversion of a 
quark-gluon plasma to hadronic matter then is usually described by the mixed phase which 
is characterized in terms of the volume fraction a, a = Vu/Vtot- Here Vh is the volume of 
the hadronic subphase in the mixed phase of volume V tot . Then any density, e.g. entropy 
density, in the mixed phase can be written as 
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s(r) = a(r)s H + (1 - a(r))sQ, 



(4) 



where sh and sq are the entropy densities in the hadronic and plasma subphases, respec- 
tively. 

If the phase transition maintains chemical equilibrium, entropy and particle densities 
can be described via thermal distributions. For a baryon free system, temperature is the 
only parameter for the thermal distributions and moreover, temperature stays constant at 
T c during the phase transition. Thus entropy densities, sh and sq, are constant during the 
phase conversion. Solving Eq. (|) for a(r) and using Eq. (||), we get the time dependence 
of a(r): 

i \ Th t ~ t Q / k \ 
a{r) = , (5) 

th -t q t 

where tq and r# are the proper times at the beginnig and end of hadronization. 

Eq. (^) implies that as the total volume of the mixed phase increases linearly in r 
according to Eq. (|3]), the plasma subvolume decreases linearly while the hadronic subvolume 
increases linearly. 

V Q = (1- a(r))V tot (r) = ^^V(r Q ) 

TH ~ Tq 

V h = a(r)V tot (r) = T JLl^^ V {r Q ) (6) 

TqTk-Tq 

Since the assumption of the chemical equilibrium during the phase conversion implies 
that the particle densities are kept constant, the number of particles hadronized should also 
increase linearly in r. This can be seen as follows: Since d A N/dx A is the particle density 
with respect to the 4-volume and dx A = rdrdyds, the density of particle species V can be 
rewritten as 

Pv(t) = 4rT" = «( r )Py ( T c) for r Q < r < r H . (7) 
rdyds 

One should note that number density in Eq. (0), py(r), is a density with respect to the 
total volume of the mixed phase, while py(T c ) is the one with respect only to the hadronic 
subvolume. Namely, pv( T ) = a ( T )Pv C^c)- From Eq. (|7|) we see that dN v /dyds indeed grows 
linearly in r. The corresponding rate equation for the hadronization of particle species V 
can be written as 

dN v 1 t h H/rn \ _ 1 



rdrdyds r th — tq t 

with 

t h 



p^(T c ) = -a v , (8) 



TH ~ Tq 



Pv(Tc), (9) 



which is constant during the mixed phase since py(T c ) is kept fixed by its value for the chem- 
ical equilibrium at a constant temperature at Tc- For a baryon-rich quark-gluon plasma, 
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dynamics will be different from Eq. @ and also py(T c ) may vary as Tc changes due to the 
reheating. 

Eq. (P) is a consequence of the gradual phase conversion but, for different hadronization 
scenarios, such as suddedn hadronization of strange particles, the result might be qualita- 
tively different. 

In the mixed phase there will also be two-body collisions among the increasing number 
of hadrons and free-decays of those hadrons. In the next section we consider them explicitly 
to set up a rate equation and solve it to obtain a time dependence of the number densities 
of a vector meson V. 



III. RATE EQUATIONS WITH SECONDARY COLLISIONS AND FREE-DECAY 

In a system of many different particle species, collisions between particles create or 
annihilate certain particle species V. V can decay freely also. Those processes change the 
number density of V on top of the hydrodynamic evolution of the system. In this section 
we consider three types of processes which changes the number of vector meson V: 

i + j^V + X, l + V^X, V^X. (10) 

These processes represent production of vector meson V through a collision between species 
i and j, annihilation of V via collisions with species /, and the free decay of V, respectively. 

The rate equation for the number density of particle species V per unit volume and time 
can be written as 

^ = — + E < ™ PiitiPii?) -E<™ >fv fi{*)pv{x) - Tfp v (x), (11) 

where the first term on the right hand side is the hadronization rate discussed in the previous 
section. The second and third terms are the gain and loss through the two body collisions 

and the last term is the loss due to the free decay as in Eq. (ITO) . In Eq. (0) the scattering 

v x — — 

cross section < crv > i '■ denotes the proper average over the momentum distribution of 

particles involved, and the density of particle species i, Pi(x), depends on time through its 
rate equation. Due to the lack of detailed knowledge, however, we will solve Eq. ( pTT|) only 
for the vector mesons of interest and assume that the number densities for other particle 
species, which appear in the right hand side of Eq. (|TTD, are determined from the chemical 
equilibrium. 

Although Eq. (|7|) looks same in each phases except for the hadronzation term, its r- 
dependence is different for the mixed and hadronic phase. Thus Eq. ( |TTD should be solved 
separately in each phases, as discussed in the following subsections. It should also be noted 
that in the mixed phase Pi(x) is the density of particle species % with respect to the total 
volume of the mixed phase, as in Eq. (|7]). 
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A. Mixed phase 



In the mixed phase, the particle densities with respect to the total volume can be written 
in terms of the densities relative to the hadronic volume as in Eq. (fil). Let us define Ay 

j \A,M 

and A v as 

\P,M sr^ „ ^V.X H H \A.M H 

Xy =2^< (TV > id Pi Pj ,X V ' =^<av > lv Pl , (12) 

where superscript M which denotes the "mixed phase" is used to distinguish similar quanti- 
ties defined for the hadronic phase. In Eq. fll2|) all the densities are taken as those relative to 
the hadronic subvolume. For the hadronization maintaing chemical equlibrium, pf = pf(T c ) 
depends only on the transition temperature, and Xy M and Xy M depend only on the collision 
geometry, but not on the time. 

Then, using the relation d 4 x = rdrdyds, we can simplfy the rate equation, Eq. (|TTD as 



Tr\^ds) =av + Ta{T)Xv - (a(T)V +Tv) ^ (13) 
Eq. ( [T3| ) can be easily solved as 

dNy , . \ A ' M t^ Ir v 4 . M iri«t\ T 

— fr) = t v T QQ~y CaA v +L v > T 

dyds 

■ [[ay + c 2 a X P y M {T ' -^ 2 ]r>-^ M r Qe (c^ M +rry dT > 

for t q <t <t h (14) 



where c a = t h /(t h - r Q ). 

In order to get insight for this complicated looking solution and get numerical values, let 
us define a ratio R v between the gain and loss terms of V due to secondary collisions and 
free-decay as 

2/ \\P,M \ p , M Ttj T ~ TQ i 

pO _ Ta ( T )V _ Ay' T H TH - TQ 1 



[a(r)Xp M +TY)t Pv (t) Xp M t Pv (t) 1 + 



V 

aX V,M 



Ry{r) is a r-dependent parameter, which controls the degree of balance between the gain 
and loss terms. For Ry{r) = 1, we expect that dN/dyds should linearly increase as a 
function of r. In our numerical estimation for the Eq. ( |14"D and thus for the number of 
dimuons, we will approximate R v {t) as its value at r = tjj. 

Ay' r H a v i _|_ v 



This approximation is mainly due to the lack of the detailed knowledge on the parameters 
involved, namely Ay M and Xy' M . However, it allows our estimation of the number of vector 
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mesons V and the corresponding dimuons produced in the mixed phase from the analysis of 
the purely hadronic experimental data. This point will be discussed later. 

Writing Ry = R v (l + par), we eliminate Ay in terms of Ry in Eq. (|Hp. 



v 



dyds 



V / s r \ A ' M T„ (r \ A ' M 4-T tot W 

-( T ) = a v r CaX v T «exp" (CQ V +r v > T 



f [1 + R v \P M c a {T '~ TQ)2 ]r'-^ M r Q eX p(-^ M +r->' dr > (17) 
Jt Q r 



B. Hadronic phase 

Worked out for the mixed phase, time evolution of particle species V in the hadronic 
phase can be obtained from Eq. (|TTD omitting the term for the hadronization, i.e. the first 
term in the right hand side. However, one should note that in the hadronic phase particle 
density is governed only from the expansion of the system and inversely proportional to r 
as in Eq. (g), while in the mixed phase linear increase in the hadronic subvolume during 
the phase conversion and constant rate of hadronization makes the particle density in the 
hadronic subvolume constant. 

From the fact that 

dN v 

- — — = T0 V {x) = const. (18) 
dyds 

in the hadronic phase and with the definition of \ H 's for the hadronic phase as 

^v' H = < av > ij X {dN / dyds) i(dN / dyds) j 

4 ,H = E<™ >fv (dN/dyds) h (19) 
i 

we obtain the rate equation for the particle density of V in the hadronic phase || as 

^(^) = hp«- { hr + nA (20) 

dr \ dyds J r t dyds 

Solution of Eq. © is 

dyds dyds Jt h t 

for T H <T < Tf, (21) 

where Ay is defined as 

Ay(r) = X V ' H In — + Tf{r - r H ). (22) 
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Similarly for the mixed phase (Eq. (|16D), we define in the hadronic phase the ratio Ry 
between gain and loss terms for particle species V due to rescattering and free-decay as 



R 



A 



P,H 



V 



(\$' H + TT v lt )(dN v /dyds) 



(23) 



Loss due to free decay of particle V is included in Eq. (|23|) , which is omitted in the definition 
of R v in Ref . . The free decay term is not negligible compared to other terms and inclusion 
of this term is consistent when we consider the effect of free decay terms. 

We will make approximation in the final calculation by taking the value for R v at r = th 
in Eq. (|23|). Then from that fact that Xy M = \y H /t^ at r = it is easy to check that 



R v defined in Eq. (|T3J) and Eq. (^) are single valued at r = r H . 

Defining Ry = R v (th)(1 + T fijr ) ; as in the mixed phase, Eq. (|2TT) becomes 

A ' 1 v 



dNy(T) dNy(T H ) 



cxp 



HM')] +RyX^ H ex^- A -^ 



dr' 



TH 



— exp 



[AvM] 



(24) 



dyds dyds 

After breakup of the system, particles decay freely and we have the exponential decay 



law 



dN^ 



dN^ 



r 



v 



dyds dyds 



(r f ) ex P [-ir (t - T f )] for r f <r. 



(25) 



C. Production of muon pairs 

Having the time dependence of the abundance of vector meson V in all phases we can 
calculate the dilepton production as 



dN l y 1 ' 

dyds 



v 



1 v 



, dN v 

dyds 

t h dN v 

dr— f (r) + 
t q dyds .i t u 



v dN v 
dr- — —it 



dyds 



1 dN v 



-ptot 
1 V 



dyds 



fa) 



(26) 
(27) 



where T l v l is the partial decay width for the the dilepton channel and Ty 1 is the total 
decay width of V. The first and second terms in Eq. ( p7|) are production in the mixed and 
hadronic phases and the last term is the number of muon pairs decayed after the freeze- 
out. It should be emphasized that in p+p collisions or when thermalized dense system is 
not formed during the heavy-ion collisions, we have only the last term in Eq. (p7|) for the 
production of muon pairs. 

Importance of the first two terms can be visualized from an idealized case when the 
gain and loss terms balance each other perfectly, namely for R v (t) = 1. In this case the 
number of particle V in the mixed phase is a linearly increasing function of t and that in 
the hadronic phase is a constant: 
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dNy . ( v T-rqdNv s , ^ ^ 

Thus we have for the number of muon pairs 



diV{, +r _ dN v . , r_^,l 
dyds dyds 



Tf (-AT M + Ar H ) :i 



(29) 



where By = T l y l /Ty 1 is the branching ratio of V in vacuum for the dilepton channel 
and Atm = (t h — tq) and At h = (t/ — t# ) are the time durations of mixed and hadronic 
phases, respectively. 

The first two terms in Eq. ( p9|) is the muon production in the mixed and hadronic phases, 
respectively, and the last term is those decayed after freeze-out. For the ideal case of perfect 
balance between gain and loss terms the production of muon pairs in the mixed and hadronic 
phases relative to those decayed after freeze-out can be simply estimated by comparing the 
magnitude of T^^Atm + At#) with 1. Please note the factor 1/2 in the contribution from 
the mixed phase which is due to the average of the linearly increasing number of vector 
mesons V. 

Having large total decay width, r* oi = 0.77 fm _1 , p mesons produce more dileptons in 
the medium than after freeze-out for relatively small lifetimes of a few fermi, while for other 
particles, e.g. for uj mesons whose total decay width is T 1 " 1 = 0.043 fm _1 the contributions 
from the mixed and hadronic phases becomes significant only when (|Atm + Ath) > 20 fm. 

For a purely hadronic system one can estimate the lifetime of the hot matter, Ath, if the 
number of vectors mesons V at freeze-out is known and the number of dileptons from them 
are measured. Or if we take the ratio of muon pair production from p and u, for which we 
expect to be 1 as in the p + p collisions due to the small mass difference and the same quark 
contents, we can get valuable information of Ath, as suggested in Ref. p. 

Even though the dependence of the muon production on the lifetime of the mixed phase 
is weaker by half than that in the hadronic phase, muon pairs decayed from the p mesons 
in the mixed phase may be the dominant one when Atq is much larger than Ar#. This 
could be the case when the hadronization transition is of first order. During the first order 
phase transition the mixed phase is expected to be long because large amount of entropy 
carried by gluons should convert to the entropy of hadrons. In this case we expect large 
enhancement of dimuons under the rho peak. 

If the mass shift of vector mesons due to restoration of symmetry occurs, those dimuons 
produced in the medium will appear at different invariant mass as suggested in ref. [ 0,0]]. 
However, what we deal in this paper is the total number of muon pairs regardless of their 
invariant mass and is not affected from the possible mass shift. 



IV. RESULTS 

Eq. (|29| ) is a result for the perfect balance between the gain and loss terms, i.e. R v {t) = 1 
for any r. For the evaluation of Eqs. (|P7D, (|24|), and (p8|) in general, values for the quantities 
Xy M , \y H and Ry(r) should be determined together with the time scales, tq, th, and r/. 
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Following the same arguement by Ref. 0, one can estimate Ay' for rho, omega and 
phi mesons from the data for J/ty and <j)/{p + to) ratio, both measured in 200 A-GeV O+U 
and S+U collisions by NA38. Underlying assumption of the arguement is that both the 
suppression of J /^ and enhancement of <p/(p + to) ratio can be simultaneouly described by 
the secondary collisions among hadrons through the rate equations, Eq. ( pip for a purely 
hadronic matter. 



When we rewrite Ay' as 



Ay H = 7vA^J, (30) 



where 



£ < crv >f Y (dNi/dyds) 
£ < <tv >u/y {dNi/dyds) 



7y involves only ratios of suitably weighted absorption cross sections and it does not depend 
on the details of the nuclear collision dynamics. From the analysis of NA38 data we have 
= 0.9 and 7y can be approximated as the ratio of cross sections involved, 7 P = 

7^ ~ 12. Thus we get X^ ,H = \^' H ~ 11. One can also estimate Ay M from the relation 
Ay M = Ay ,if /rff, as both quantities are approximately independent of r. 

As we have defined the quantity Rv{r) in the previous section we can eliminate Ay' M , 
and Xy H in terms of i?y(r). We further approximate R v (t), which is a function of r, with 
its value at at r = th, i-e. Ryij) ~ R v {jh) and take it as a parameter which controlls the 
degree of equilibration. Namely for Ry > 1 creation due to secondary collisions wins over 
the annihilation by the secondary collisions and free decay, and for R v < 1, vice versa. 

Time scales involved, tq, th and Tf, should be determined from the hydrodynamic 
equations together with a model equations of state for each phase. Instead of making 
models we rather satisfy ourselves using those time scales as variables. 

It is convenient to define a ratio 

= dN v (r)/dyds 
V dN v {r H )/dyds' 1 ' 

In Fig. 1 we have plotted Yy for p and uj mesons with arbitrarily chosen values, tq = 5 
fm, At m = (th — tq) = 20 fm, Ar H = (r/ — t h ) = 35 fm. Solid lines are the results of 
the approximation R v (r) ~ R v (t h ) = 1 while the dotted lines are for the ideal case, i.e. 
Eq. fl2"9"|). The approximation Ryij) ~ R v {jh) = 1 underestimates the loss due to free decay 
and one should note the importance of this contribution especially for the rho mesons. 

The unknown overall factor [cINv(th) / dyds] can be removed by taking a ratio of muon 
production from p to u mesons. In p + p collisions p and uj mesons are produced in equal 
numbers, which can be understood from the small mass difference and same quark contents. 
Execept for the medium effect on the number of the two vector mesons, it seems natural to 
approximate that [dNyiju) / dyds] for p and to are same. Then we get from Eq. (|29|) for the 
perfect balance 

x = dN l p +l -/dyds = B p [t^ 2 At m + At h ) + i] 



dN l p +l ~ 


/dyds 


B P 


r*°*(|Ar M + At h ) + 1 


dm +l - 

UJ 


' /dyds 




T™{\At m + At h ) + 1 
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with B p /B w = 0.65. 

In Fig. 2 we show calculated X as a function of the lifetime of the mixed phase, Atm 
(t h — tq) for two fixed values of Ar H = (t/ — r H ). 

When At m /2 + Ar H <C 20 fm, we can approximate X as 



(34) 



which is a linear function of Atm/2 + Ar#. 

Taking At m = in Eqs. (|33| ) we get the ratio, X for a purely hadronic system. 



Bp 

X 



T^Ath + 1 



\V^At h + 1] 



0.65(rf Arff + 1). (35) 



According to Eq. (|35|) measurement of X larger than 0.65 implies the formation of a hot 
and dense medium in thermal equilibrium. If the hadronization of a quark-gluon plasma 
occurs as a first order phase transition, value of X much larger than 0.65 is expected. 

For a purely hadronic system we get from Eq. (|35|) Ath by measuring the ratio of 
dileptons from p to u mesons. Thus careful anaysis of muon production in the p and to 
channels may reveal valuable informations even on the formation of the hadronic fireball as 
suggested by Heinz and Lee f|. In Ref. || the dependence of X on the total transverse 
energy, Et,is investigated. In this paper, however, we merely note the fact that as the beam 
energy or the total transverse energy increases, the initial temperature of the hadronic fireball 
increases and so does Ar#. 

For the beam energy just above that needed to make a deconfmement transition, a state 
of mixed phase with small Atm will be formed and as the beam energy increases, so rapidly 
does At m with more or less saturated Ar H . The dependence of (Atm/2 + Ar H ) on the 
beam energy will be different from that of At# at lower energies. 

For high enough beam energies a quark-gluon plasma will be made and Atm and Ath 
should be functions of the initial temperature of the QGP. Thus as the beam energy increases 
there will be two characteristic changes in the behavior of X at the energies of barely forming 
a state of mixed phase and a quark-gluon plasma. 

As in the case of the average hadronic spectra as a function of the beam energy, 
characteristic changes in slopes of X as a function of the beam energy can serve as a signature 
of QGP. Further studies are needed to investigate how drastic will be the change. 

Under the debate is the change of mass of vector mesons in the hot and dense medium 
due to the restoration of chiral symmetry, as discussed in the last section. If the mass of 
vector mesons decreases due to chiral retoration in the hot and dense medium, dileptons 
decayed during the mixed phase will peak at smaller invariant mass. 

In this case we may have in the invariant mass spectrum a new peak from the mixed 
phase at the mass correspoding to the hadronization transition, my(Tc, Pb) together with a 
peak due to the free decay after freeze-out. Dileptons decayed in the hadronic phase whose 
temperature and baryon density change continuously will spread out in the mass region 
between mv(T c , Pb) and m v in the free space and can hardly be distinguished from the 
background. However, it should be stressed that even when the mass of vector mesons V 
changes in the hot and dense medium, our results so far are valid as we have considered the 
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total number of muons pairs decayed from vector mesons V in each stages of the system but 
not the invariant mass spectrum. 

There may also a change of the decay width of vector mesons in the hot and dense 
medium [jSHltiH. In this case we need to put r- dependent Y^ij) in every equations and the 



relation between r and the temperature and baryon density should be determined from the 
dynamics of the fireball. 

V. SUMMARY 

We explore interesting aspects of muon pair production in the relativistic heavy-ion 
collisions and have shown that careful measurement of muon pairs under p and o> peaks may 
reveal valuable informations on the formation of quark-gluon plasma. 

Especially, the ratio of muon pairs from p to u mesons, X, carries the information on 
the lifetime of a fireball either in the hadronic or quark-gluon plasma phase. Characteris- 
tic changes in X as a function of increasing beam energies should serve as a signature of 
formation of a quark- gluon plasma. 

Analysis of muon pair production decayed from vector mesons is made by solving rate 
equations which is constructed by considering a hadronization rate, creation and annihilation 
of particles via secondary collisions, and free-decays. For the hadronization rate of vector 
mesons of interest we consider a baryon-free quark-gluon plasma which is invariant under 
Lorentz boost, and assume equilibrium hadronization at a contant temperature. This implies 
a constant rate of hadronization during the mixed phase. 

However, inside a hot medium broadening of mass due to frequent collisions will spread 
the muon peaks under rho mesons and this broadening of rho peak in the dimuon spec- 
trum may offset the increase of the rho peak to make to measurement of the increase hard 
Studies of the collisional broadening are needed for any realistic estimation of the 



enhancement of dimuons under the rho meson peak from the medium. 
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FIGURES 



Fig.l: Number of vector mesons per unit rapidity and transverse fucntion of r. 

Solid lines are the results of the approximation R v (r) ~ R v {jh) = 1 while the dotted lines 
are for the ideal case, i.e. Eq. (|29|). 



Fig. 2: Ratio of muon production from p to w mesons as a function of lifetime of the 
mixed phase. Here time length of hadronic phase is fixed as 5 fm. 
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